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Abstract 

In this work, we calculate the rare decays D*" and B~^ — > D*'^T^ 

in perturbative QCD approach with Sudakov resummation. We give the branching 
ratio of 10~^ for B^ D*^ , which will be tested soon in B factories. The decay 
B+ D*+K has a very small branching ratio at 0(10 due to the suppression 
from CKM matrix elements It may be sensitive to new physics contributions. 
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1 Introduction 



Perturbative QCD (PQCD) method for B decays has been developed for some years pQ. It 
is successfully applied to exclusive B meson decays recently, such as B ^ nn [2,, i? — > Kn 
jS], i? vrp [?j and other channels |5]. 

Very recently, pure annihilation type B decays are discussed in the PQCD approach, such 
as 5° D-K+ and B+ D+K^ 013 decays and B+ D+(f) decay [8 . It is found that 
the annihilation type decay B^ DjK^ has a sizable branching ratio of 10~^, which has 
already been measured by experiments 0. In this paper, we will continue to compute the 
branching ratios of similar decays 5° — ^ D*~K~^ and B^ D*^K^ . 

Because the four valence quarks in the final states D* and K{K) are different from the 
ones in the B meson, the rare decays B ^ D*K are pure annihilation type decays. In the 
usual factorization approach (FA) [Hj, this decay picture is described as B meson annihilating 
into vacuum and D* and K meson produced from vacuum afterwards. To calculate this decay 
in the FA, one needs the — > K form factor at very large timelike momentum transfer 
0{Mb)- However the form factor at such a large momentum transfer is not known in FA. 
The annihilation amplitude is a phenomenological parameter in QCD factorization approach 
(QCDF) PU], and the QCDF calculation of these decays is also unreliable. 

In this paper, we will calculate these decays in PQCD approach. Similar to the B DsK 
decays, the W boson exchange induce the four quark operator hd —>■ cu or hu — > dc, and 
the ss quarks included in D^K are produced from a gluon. This gluon attaches to any one 
of the quarks participating in the four quark operator. In the rest frame of B meson, the 
produced s or s quark included in D*K final states has (9(Mb/4) momentum, and the gluon 
producing them has ~ (9(M^/4). This is a hard gluon, so we can perturbatively treat 
the process where the four-quark operator exchanges a hard gluon with ss quark pair. It is 
a perturbative six quark interaction now. 

In PQCD, the decay amplitude is separated into soft ($), hard (iJ), and harder (C) 
dynamics characterized by different scales. It is conceptually written as the convolution. 

Amplitude- jd^k^d^k^d^ i:i[C{t)^B{ki)^Dt{k2)^K{h)H{ki,k2,k^,t)\, (1) 

where fcj's are momenta of light quarks included in each meson. C{t) is Wilson coefficient 
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which results from the radiative corrections at short distance. $m is the wave function which 
describes hadronization of the quark and anti-quark into the meson M. H describes the four 
quark operator and the quark pair from the sea connected by a hard gluon whose scale is at 
the order of Mb, so the hard part H can be perturbatively calculated. 

Some analytic formulas for the decay amplitudes of 5 — > D^K decays will be given in 
the next section. In section El we give the numerical results and discussion. Finally, we 
conclude this study in section El 



2 Analytic formulas 

We consider the B meson at rest for simplicity. In the light-cone coordinate, the B meson 
momentum Pi, the meson momentum P2 and K meson momentum P3 are taken to be: 

Mb, X Mb, o , Mb, 9 , 

Pi = 1, Ot), P2 = ^(1, r^ Ot), P3 = ^(0, 1 - t\ Ot), (2) 

where r = M£,*/Mb and we neglect the K meson's mass Mk- The meson's longitudinal 
polarization vector e is given by e = -7^^ — {'^i~'^'^i^t)- Denoting the light (anti-)quark 

^ s 

momenta in P, D* and K mesons as fci, A;2, and fcs, respectively, we can choose ki = 
(xiPi^, 0, ki-r), k2 = {x2P^ ,0,'k2T), ^3 = (0, xsPg", ks^). The decay amplitude in eq.((H) 
leads to: 



Amplitude ~ J dxidx2dx3bidbib2db2bsdbs 

Tr [C{t)<l>B{xu 6i)<l>D| {X2, 62, e)<l>K{x3, b^)H{x.,, b„ e, t)St{x,) e^^W] , (3) 

where bi is the conjugate space coordinate of kix, and t is the largest energy scale in H. The 
large logarithms \n{mw/t) coming from QCD radiative corrections to four quark operators 
are included in the Wilson coefficients C{t). The last term, e"'^'-*-', contains two kinds of 
logarithms. One of the large logarithms is due to the renormalization of ultra-violet diver- 
gence Intfe, the other is double logarithm ln^6 from the overlap of coUinear and soft gluon 
corrections. This Sudakov form factor suppresses the soft dynamics effectively 11 . Thus it 
makes perturbative calculation of the hard part H applicable at intermediate scale, i.e.M^ 
scale. 
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The B meson wave function for incoming state and the D* and K meson wave function 
for outgoing state with up to twist-3 terms are written as: 

^bM^, b) = [( A75)a/3 + MB75a/3] M^, b) , (4) 



<fDj,a/3(x,6) = [/(A + Md.)„;3] ^Dtix,b), (5) 



+ "^oif75(/' - l)0]^(a;3,&3) , (6) 

where f = (0, 1, Ot) oc P3, n = (1, 0, Ot), and = Mj^/irriu + ms). 

Since the wave functions are process independent, one can use the same forms constraint 
by other decay channels [5] to make predictions here. Now the hard part H in eq.® is the 
only channel dependent part for us to calculate perturbatively. 

2.1 Dl K+ decay 

In the decay D*~'K^ , the effective Hamiltonian at the scale lower than Mw is given 

as [ig 

i/eff = ^KtKd [CA^^)0,{^^) + CM02{^^)] , (7) 
Ol = {bd)Y^A{uc)v-A, O2 = {bc)v~A{ud)Y^A, (8) 

where Ci^2{fJ') are Wilson coefficients at renormalization scale /z, and summation in SU(3)c 
color's index a and chiral projection, X]a^a7^(l ~ 75)5'q: abbreviated to {qq')v~A- The 
lowest order diagrams of B^ D*^K^ are drawn in FigQ according to this effective Hamil- 
tonian. There are altogether four quark diagrams. Just as what we said above, the decay 
B D*K only has annihilation diagrams. 

By calculating the hard part H at the first order of a.^, we get the following analytic 
formulas. With the meson wave functions, the amplitude for the factorizable annihilation 
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(d) 



Figure 1 



Figure 1: Diagrams for —>■ D* decay. The factorizable diagrams (a),(b) contribute 
to f!^'\ and nonfactorizable (c), (d) do to Ma. 



diagrams in Fig^^a) and (b) results in Fa as: 



Jo Jo 

X { (1 - r^) (1 - (1 - r^)x3) 0^(X3, 63) + r (1 - r2)x3r^0^(a;3, ^3) 
- r(l + r2)r^0^(x3, h)}E}{tl)ha{x2, X3, 62, ^3) - { (1 - r2)x20^(x3, h) 

+ 2r(l - - a;2)r^0^(x3, 63) }^}(ta)/^a(l - ^3, 1 - X2, 63, ^2) 



(9) 



In the function, Cp = 4/3 is the group factor of SU(3)c, and = rnQK/Ms- The functions 
E'jr, tl^"^, ha will be given in the appendix. The distribution amplitudes (pM, are given in the 
next section. The amplitude for the nonfactorizable annihilation diagrams in FigGJc) and 
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(d) results in 

2 rl poo 

Ma = —^===64:7:0 fM'^ / dxidx2dx^ / 62'^&2 &i)0D|(a;2, ^2) 

X {(l-r^) ((l-r2)(l-x3)-r\2)0:^(x3,&2) 
- r (x2 - (1 - r^) (1 - X3)) rK(l)Ki^3, h) 
- r (x2 + (1 - r^)(l - X3)) rx0x(3;3,&2)}-Em(t;tt)^aH3;i,a;2,a;3,6i,62) 
-{(1-r^) {{1 - r')x2 + r') <Pi{xs,h) 
+ r{x2 - (1 - - X3))r;^0]^(x3, 62) 
+ r{x2 - 2r2 - (1 - r2)(l + X3))r Kcf^x^, b2)}Emitl)h^^\xu X2, X3, 61, 62)] , (10) 

where xi dependence in the numerators of the hard part are neglected by the assumption 

Xi < X2, X3. 

Comparing with D~K^ [6^, we find that the leading twist contribution, which is 

proportional to 0;^, is almost the same. However the sub leading twist contribution, which 
is proportional to r r^, change significantly, especially for the 0^(x3, 62) terms in (jH)) and 

The total decay amplitude for D*~K^ decay is given as A = JeF!,''^ + Ma. The 

decay width is then 

(^2 A^S 

T{B' = ^1^(1 - r')\V:,V^,A\'. (11) 

The charge conjugate decay 5^ D*~^K~ is the same value as 5° D*~ just because 
iKfe^dl is the same as iKfeV^dl- Since there is only one kind of CKM phase involved in the 
decay, there is no CP violation in the standard model for this decay channel. 

2.2 ^ Dl+lf decay 

The effective Hamiltonian related to B^ D*^l{^ decay is given as ^21 

i/eff = ^Kft^cd [Cl{fi)0,{fi) + C2(/i)02(/i)] , (12) 

Oi = (bd)v~A{cu)v-A, O2 = (bu)v-A{cd)v~A- (13) 
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The amplitude for the factorizable annihilation diagrams results in —Fa ^\ The amplitude 
for the nonfactorizable annihilation diagrams is given as 

M'^ = -^==M%CfM% I dxidx2dx3 / hidhih2dh2 6i)0d|(x2, 62) 

V ^J^c Jo Jo 

X {(1 - 2r^)x2(^i{x3, 62) + r (x2 - (1 - r^){l - X3)) rK4>^{x3, 62) 
+ r (x2 + (1 - r2)(l - X3)) rK<p]i{x3, h2)]E'^{t]^)h^a\xi^ ^2, X3, 61, 62) 
- {(1 - r^) (1 - X2r2 - (1 - r2)x3) 0^(x3, &2) 
- r (x2 - (1 - r^)(l - X3)) rx0x(3^3, ^2) 
-r{x2- 2r2 - (1 - r2)(l - xg)) r^0^(x3, 62)}^;„(4)/if (^^i, 2:2, X3, 61, 62)! • (14) 



Comparing with the decay to two pseudo-sealer final states —> D^K p|, the change 
only occur in terms. Therefore, the non-factorizable contribution in decay does not 
change much. 

The total decay amplitude A' and decay width F for B^ D*^K^ decay are given as 

A' = -fnFj,'^ + (15) 

r{B^ Dl+lf ) = ^^(1 - r^)\V:,V^,Af. (16) 

The decay width for CP conjugated mode, B^ D*~K^, is the same value as B~^ D*^K^ . 
Similar to the B^ decay, there is also no CP violation in this decay within standard model. 



3 Numerical evaluation 



We use the wave function of B meson written as [S] 

0b(x, b) = Nbx'^{1 — xY exp 



M| x^ 1 



(17) 



where Nb is a normalization factor. Because the mass difference between and Z)* is not 
large, for simplicity, their wave functions are chosen to be the same OE]: 



0£,*(x,6) 



:/d|x(1 - x){l + aD|(l - 2x)}exp 



Since c quark is much heavier than s quark, this function is peaked at c quark side, i.e. small 
X region. 



The K meson wave functions are given as 

^i{x) = ^^6x(l-x)|l-af -Se + af -^(-l + Sa}, (19) 

rAx) = ^;^{i+<i-^(-i+3a+<2-^(3-3oe^+35a}, (20) 

<f>l{x) = ^-^(l-2a;){l + a^-3(-3 + 5a}, (21) 

where ^ = 2x — 1. The parameters of these distribution amphtudes calculated from QCD 
sum rule [121 given as 

af = 0.17, = 0.2, af^ = 0.212, a^^ = -0.148, = 0.0527, (22) 

for moK = 1.6 GeV. In addition, we use the following input parameters [U]: 

fe = 190 MeV, Jk = 160 MeV, fot = 241 MeV, (23) 
moK = 1.6 GeV, tOb = 0.4 GeV, a^j = 0.3, ud* = 0.3 GeV. (24) 

For branching ratio estimation, we use the CKM matrix elements and the lifetimes of B 
mesons as following 

\Vud\ = 0.9734, \Vub\ = 3.6 x 10-^ (25) 

\Vcb\ = 41.2 X 10-^ \Vcd\ = 0.224, (26) 

tb± = 1.67 X 10"^^ s, tbo = 1.54 x 10"^^ s. (27) 

The predicted branching ratios are 

Br(fi° ^ D*-K+) = 2.7 x 10'^ (28) 
Br{B+ Dl+lf ) = 4.0 x 10"^ (29) 

The branching ratio of D*^l{^ is much smaller than that of 5° — > D*~K~^, due to 

the suppression from CKM matrix element | V*fycd \« iKbKtdl- C-H Chen has also given 
the branching ratio of D*^K^ [7\, our result agrees with his. 

In these processes of i? ^ D*K, only the longitudinal polarization of D* has contribution. 
If we ignore the difference between Jd^ and fo* and the difference between M^s and M^*, 
the branch ratios are thought to be the same as the corresponding B DgK decays. But 

8 



in fact the branching ratio of D*^K~^ is a httle smaller than that of DjK~^, 

because the contribution of twist-3 wave function in B^ D*~K^ decay become negative. 
Although the larger makes the branching ratio larger, a larger M^). leads to a smaller 
branching ratio. The effect of Md* is more dominant than that of fo* ■ In decays of charged 
B"^, fn* makes branching ratios larger, and M^* has little effect, such that the B^ D*^K^ 
decay in this case has a little larger branching ratio than B^ D^K^. 

For the experimental side, there are only upper limits given at 90% confidence level ^3] , 

Br(E° ^ D*-K+) < 1.7 X 10-\ (30) 
Br(5+ ^ D*+K^) < 1.1 X 10"^ (31) 

Obviously, our results are consistent with the data. 

In addition to the perturbative annihilation contributions, there is also a hadronic picture 
for the B^ D*^K^. The B meson decays into D~ and p^, the secondary particles then 
exchanging a K^, scatter into D*~ , through final state interaction afterwards. For 
charged decay, B^ meson decays into D^p^ then scatter into D*^ and by exchanging 
a K~^. But this picture can not be calculated accurately. In [6^, the results from PQCD 
approach for B^ D^K^ decay were consistent with the experiment, which shows that the 
soft final state interaction may not be important. 

The calculated branching ratios in PQCD are sensitive to various parameters, such as 
parameters in eqs. ()22ll!Mj) . It is necessary to give the sensitivity of the branching ratios 
when we choose the parameters to some extent. Tabled shows the sensitivity of the branch 
ratios to 30% change of an*, ud* and Uf,. Here we don't present the sensitivity to a^p2T 
because the branching ratios are insensitive to them. It is found that the uncertainty of 
the predictions on PQCD is mainly due to ujd*, which characterizes the shape of wave 
function. 

From the above discussions, we can derive the uncertainties of the branching ratios within 
the suitable ranges on ud*, an* and Ub- The branching ratios normalized by the decay 
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Table 1: The sensitivity of 



the decay branching ratios to 30% change oi cud*, cld* and cub. 
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[W ) 
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Br(S+ ^ D*+lf) 


0.21 




3.17 


4.66 


0.30 




2.65 


3.95 


0.39 




2.11 


3.20 




Br(S° 




Br(S+ ^ 


0.21 




2.40 


3.68 


0.30 




2.65 


3.95 


0.39 




2.92 


4.24 




Br(S° 




Br(S+ ^ £>:+^°) 


0.21 




3.12 


3.87 


0.40 




2.65 


3.95 


0.52 




2.31 


4.18 
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constants and the CKM matrix elements result in 



Br(S'^ ^ D*-K+) = (2.7 ± 0.6) x IQ-' ( ^^^^ 1 ( ^Y^l}^ ^ (32) 

^ ' ' ^ > Vl90MeV-240 MeVy VO-0412 ■ 0.9734; ' ^ ' 

Br(B+ ^ L>:+X°) = (4.0 ± 0.8) x 10'^ ( ^ ' ( -JYmlM^] ' . (33) 

^ s ^ V J \^190 MeV-240 MeVy VO-0036 • 0.224/ ^ ^ 

4 Conclusion 

In this paper, we calculate the —>■ D*~K^ and —>■ decays in PQCD approach. 

These two decays occur purely via annihilation type diagrams because the four quarks in final 
states are not the same as the ones in B meson. We argue that the soft final state interaction 
may be small in B decays. The PQCD calculation of annihilation decays is reliable. The 
branching ratio for B^ D*~K^ decay is sizable at order of 10^^, which can be measured 
in the current B factories Belle and BABAR. The small branching ratio of B^ —>■ D^'^K^ 
(10~^) predicted in the SM, makes it sensitive to new physics contributions, which may be 
studied in the future LHC-B experiment. 
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A Some formulas used in the text 

The function E^jr, E^n, and E'^ including Wilson coefficients are defined as 

£;}(i) = a,(t)a,(t)e-^-W-^-W, (34) 

Emit) = C2{t)a,{t) e-SBit)-Snit)-s.(t)^ (35) 

E'^{t) = Ci{t)as{t) e-5«W-^^(*)-^^(*), (36) 

where 

a,{t) = C\{t)/3 + C^it), a2{t) = C,{t) + C2(t)/3, (37) 
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and Sb, Sd, and Sk result from summing both double logarithms caused by soft gluon 
corrections and single ones due to the renormalization of ultra-violet divergence. Those 
factors are given in ref . jSl El • In the numerical analysis we use leading logarithms expressions 
for Wilson coefficients Ci,2 presented in ref. |12[ 1^. 

The functions ha, h^a \ and h^a^ in the decay amplitudes consist of two parts: one is the 
jet function St{xi) derived by the threshold resummation [T3], the other is the propagator of 
virtual quark and gluon. They are defined by 

K{x2,X3,b2,b,)=St{l-x,) (^y^ Hl,'\MBy/{l-r^)x2{l-x,)b2) 

X |i7i^)(MB V(l - r2)(l - X3) b2)MMBVa - - ^3) bs)9ib2 - b,) + (62 ^ 63)} , 

(38) 

hl^\xi,X2,X3,bi,b2) = 

|^H;,^)(Mbv/(1 - r2)x2(l - X3) 6i)Jo(Mbv/(1 - r^)x2{l - X3) &2)^(&i - 62) 

Ko(MbF(,)6i), for F2.)>0 
fH«(MB,/f^6i), for <0^ 



+ {b,^b2)\x\ . , (39) 

where Hq^^(z) = Jo(-2) + i Yo(-2), and -F(j)S are defined by 

Ffi) = (1 - r2)(xi - X2)(l - xs), = xi + X2 + (1 - r2)(l - xi - X2)(l - X3). (40) 

We adopt the parametrization for St{x) of the factorizable contributions, 

^,(x) = ^^i/-\i[x(l-x)r, c = 0.3, (41) 
VttI (1 + cj 

which is proposed in ref. |16J. The hard scale t's in the amplitudes are taken as the largest 
energy scale in the hard part H to kill the large logarithmic radiative corrections: 



tl = max(MBv/(l-r2)(l-a;3), 1/^2, 1/&3), (42) 
tl = max(MsV(l-^')a;2, 1/&2, 1/&3), (43) 



C = max(MB Ms v/(l - ^')a;2(l - X3), I/61, 1/62). (44) 
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